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, Abstract. Cosmological observations have in the past few years become an increas- 

■ ingly powerful method for determining parameters within the neutrino sector, such as 

^SJ ' the presence of sterile states and the mass of neutrinos. I review the current status of 

, the field in light of recent measurements of the cosmic microwave background by the 

O ' WMAP collaboration, as well as current large scale galaxy surveys. 

O: 

^ ■ 1 Introduction 

^ ' Neutrinos exist in equilibrium with the electromagnetic plasma in the early uni- 

, verse, until a temperature of a few MeV. At this point the weak interactions 

' freeze out and neutrinos decouple from the plasma. Shortly after this event, 

, the temperature of the plasma falls below the electron mass, and electrons and 

' positrons annihilate, dumping their entropy into the photon gas. This heats 

. the photon gas while having no effect on neutrinos, and the end result is that 

' the photon temperature is larger than the neutrino temperature by the factor 

, (ll/4)^/'^ ~ 1.40. Since the present day photon temperature has been measured 

pi 1' with great accuracy to be 2.728 K, the neutrino temperature is known to be 

Q , 1.95 K, or about 2 x 10"'' eV. Since the heaviest neutrino has a mass of at least 

' about 0.04 eV it must at present be extremely non-relativistic and therefore acts 

, as dark matter. The contribution of a single neutrino species of mass m^, to the 

' present day matter density can be written as |1I2I3| 

: - 

' so that for a neutrino mass of about 30 eV, neutrinos will make up all of the 
dark matter. However, this would have disastrous consequences for structure 
formation in the universe, because neutrinos of eV mass have very large free 
streaming lengths and would erase structure in the neutrino density on scales 
smaller than Ifs — 1 'm~ly Gpc completely. This leads to an overall suppression 



of matter fluctuations at small scales, an effect which is potentially observable. 
1.1 Absolute value of neutrino masses 



The absolute value of neutrino masses are very difficult to measure experimen- 
tally. On the other hand, mass differences between neutrino mass eigenstates, 
(mi, TO2, can be measured in neutrino oscillation experiments. Observations 
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of atmospheric neutrinos suggest a squared mass difference of 6m^ ~ 3 x 10^^ 
eV^ | 4I5I6| . While there are still several viable solutions to the solar neutrino 
problem from solar neutrino observations alone, the large mixing angle (LMA) 
solution gives by far the best fit with Sm^ ~ 5 x 10^^ eV^ |7l8j . Recently the 
KamLAND reactor neutrino experiment has announced a positive detection of 
neutrino oscillations indicating that the LMA solution is indeed correct 

In the simplest case where neutrino masses are hierarchical these results 
suggest that toi ~ 0, m2 ~ <5TOsoiar, and ma ^ (5matmosphcric- If the hierarchy 
is inverted |1UI11I12I13I14I15) one instead finds ms ^ 0, TO2 ~ ^m-atmosphoric, 
and mi ~ <5TOatmospheric- Howcver, it is also possible that neutrino masses are 
degenerate |ieill7ll8ll9l2()l21l22l23l24<^^ . mi - TO2 - » (5matmospheric, 
in which case oscillation experiments are not useful for determining the absolute 
mass scale. 

Experiments which rely on kinematical effects of the neutrino mass offer the 
strongest probe of this overall mass scale. Tritium decay measurements have 
been able to put an upper limit on the electron neutrino mass of 2.2 eV (95% 
conf.) 122). However, cosmology at present yields an even stronger limit which is 
also based on the kinematics of neutrino mass. As discussed before any structure 
in the neutrino density below the free-streaming scale is erased and therefore the 
presence of a non-zero neutrino mass suppresses the matter power spectrum at 
small scales relative to large scale, roughly by AP/P ~ — 8i7y//2m |28'. 

This power spectrum suppression allows for a determination of the neutrino 
mass from measurements of the matter power spectrum on large scales, as well 
as the spectrum of CMB fluctuations. This matter spectrum is related to the 
galaxy correlation spectrum measured in large scale structure (LSS) surveys via 
the bias parameter, &^(fc) = Pg{k)/Pm{k). Such analyses have been performed 
several times before |29I3(J) . most recently using data from the 2dFGRS galaxy 
survey |31l32l33j . These investigations found mass limits of 1.5-3 eV, depending 
on assumptions about the cosmological parameter space. 

In a seminal paper it was calculated by Eisenstein, Hu and Tegmark that 
future CMB and LSS experiments could push the bound on the sum of neutrino 
masses down to about 0.3 eV |2H|- The prospects for an absolute neutrino mass 
determination was discussed in further detail in Ref. 39^ where it was found 
that in fact the upper bound could be pushed to 0.12 eV (95% conf.) using data 
from the Sloan Digital Sky Survey and the upcoming Planck satellite. 

More recently the new WMAP data, in conjunction with large scale structure 
data from 2dFGRS has been used to put an upper bound on the sum of all 
neutrino species of J2 — 0-^ (95% conf.) pi] . 

However, the exact value of this upper bound depends strongly on priors on 
other cosmological parameters, mainly Hq. In the present paper we calculate 
the upper bound on ^ rrii, from present cosmological data, with an emphasis on 
studying how the bound depends on the data set chosen. 

In addition to their contribution to the cosmological mass density neutrinos 
also contribute to the cosmological energy density around the epoch of recombi- 
nation. Neutrinos which have mass smaller than roughly 3Ti.oc, where Tree — 0.3 
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eV is the temperature of recombination, will act as fully rclativistic particles 
when it comes to CMB and large scale structure. 

In the standard model there are three light neutrino species with this prop- 
erty. However, these particles are not necessarily in an equilibrium Fermi-Dirac 
distribution with zero chemical potential. It is well known that the universe 
contains a non-zero baryon asymmetry of the order 77 — "^^ ^ 10""'^°. A 
neutrino asymmetry of similar magnitude would have no impact on cosmology 
during CMB and LSS formation, but since the neutrino asymmetry is not di- 
rectly observable it could in principle be much larger than the baryon asymmetry. 
Such a neutrino asymmetry would effectively show up as extra rclativistic energy 
density in the CMB and LSS power spectra. 

Another possibility for extra relativistic energy is that there are additional 
light species beyond the standard model which have decoupled early (such as 
the graviton or the gravitino). 

From the perspective of late time evolution at T < 1 MeV it is customary to 
parametrize any such additional energy density in terms of Ny the equiv- 
alent number of neutrino species. In Section III we discuss bounds on from 
the present WMAP and 2dFGRS data, combined with additional information 
on other cosmological parameters from the Hubble HST key project and the 
Supernova Cosmology Project. 

However, as will be discussed later, a non-zero neutrino chemical potential 
can have an effect on big bang nucleosynthesis which is profoundly different from 
simple relativistic energy density if it is located in the electron neutrino sector. 

Another important point is that any entropy production which takes place 
after BBN, but prior to CMB formation will only be detectable via CMB and 
LSS observations. One such example is the decay of a hypothetical long-lived 
massive particle at temperatures below roughly 0.01 MeV. 

2 Likelihood analysis and data sets 

The extraction of cosmological parameters from cosmological data is a difficult 
process since for both CMB and LSS the power spectra depend on a plethora 
of different parameters. Furthermore, since the CMB and matter power spectra 
depend on many different parameters one might worry that an analysis which is 
too restricted in parameter space could give spuriously strong limits on a given 
parameter. 

The most recent cosmological data is in excellent agreement with a flat 
ylCDM model, the only non-standard feature being the apparently very high 
optical depth to reionization. Therefore the natural benchmark against which 
non-standard neutrino physics can be tested is a model with the following free 
parameters: flra^ the matter density, the curvature parameter, the baryon 
density, i?Oj the Hubble parameter, n^, the scalar spectral index of the primordial 
fluctuation spectrum, r, the optical depth to reionization, Q, the normalization 
of the CMB power spectrum, fe, the bias parameter, and finally the two param- 
eters related to neutrino physics, and Ni,. The analysis can be restricted 
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to geometrically flat models, i.e. fl = flm + f^A = 1- For the purpose of actual 
power spectrum calculations, the CMBFAST package can be used. 

2.1 LSS data 

At present, by far the largest survey available is the 2dFGRS ST^ of which 
about 147,000 galaxies have so far been analyzed. Tegmark, Hamilton and Xu 
|3H1 have calculated a power spectrum, P{k), from this data, which we use in 
the present work. The 2dFGRS data extends to very small scales where there 
are large effects of non-linearity. Since we only calculate linear power spectra, 
we use (in accordance with standard procedure) only data on scales larger than 
k = 0.2h Mpc~^, where effects of non-linearity should be minimal Making 
this cut reduces the number of power spectrum data points to 18. 

2.2 CMB data 

The CMB temperature fluctuations are conveniently described in terms of the 
spherical harmonics power spectrum C; = (|aim|^), where ^^{6,<j)) — X]/m '^'m-^" 
Since Thomson scattering polarizes light there are additional powerspectra com- 
ing from the polarization anisotropics. The polarization can be divided into a 
curl-free (E) and a curl (B) component, yielding four independent power spectra: 
Ct,i,Ce,i,Cb,i and the temperature .E-polarization cross-correlation Cte.i- 

The WMAP experiment have reported data only on Ct,i and Cte.i, as de- 
scribed in Ref. 'WWWW^ 

We have performed the likelihood analysis using the prescription given by 
the WMAP collaboration which includes the correlation between different Cj's 
| 40I34I41> 42 43 . Foreground contamination has already been subtracted from 
their published data. 

In parts of the data analysis we also add other CMB data from the compila- 
tion by Wang et al. which includes data at high /. Altogether this data set 
has 28 data points. 

3 Numerical results 
3.1 Neutrino masses 

The analysis presented here was originally published in Ref. 145' , and more details 
can be found there. 

We have calculated as a function of neutrino mass while marginalizing 
over all other cosmological parameters. This has been done using the data sets 
described above. In the first case we have calculated the constraint using the 
WMAP Ct.i combined with the 2dFGRS data, and in the second case we have 
added the polarization measurement from WMAP. Finally we have added the 
additional constraint from the HST key project and the Supernova Cosmology 
Project. It should also be noted that when constraining the neutrino mass it has 
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in all cases been assumed that N^, is equal to the standard model value of 3.04. 
Later we relax this condition in order to study the LSND bound. 

The result is shown in Fig. 1. As can be seen from the figure the 95% con- 
fidence upper limit on the sum of neutrino masses is ^rrii, < 1.01 eV (95% 
conf.) using the case with priors. This value is completely consistent with what 
is found in Ref. 021 where simple Gaussian priors from WMAP were added to 
the 2dFGRS data analysis. For the three cases studied we find the following 
limits: 

Y,m^< 1.01 eV for WMAP+Wang+2dFGRS+HST+SNTa 

X; < 1.20 eV for WMAP+Wang+2dFGRS 

Y,m^< 2-12 eV for WMAP+2dFGRS 

In the middle panel of Fig. 1 we show the best fit value of i?o for a given 
f2^h'^. It is clear that an increasing value of ^rrii, can be compensated by a 
decrease in Hq. Even though the data yields a strong constraint on f2mh^ there 
is no independent constraint on ilm in itself. Therefore, an decreasing Hq leads 
to an increasing Sim- This can be seen in the bottom panel of Fig. 1. 

When the HST prior on Hq is relaxed a higher value of ^ is allowed, in 
the case with only WMAP and 2dFGRS data the upper bound is il^h'^ < 0.023 
(95% conf.), corresponding to a neutrino mass of 0.71 eV for each of the three 
neutrinos. 

This effect was also found by Elgar0y and Lahav "IB" in their analysis of the 
effects of priors on the determination of ^ m,^. 

However, as can also be seen from the figure, the addition of high-^ CMB 
data from the Want et al. compilation also shrinks the allowed range of ^ 
significantly. The reason is that there is a significant overlap of the scales probed 
by high-^ CMB experiments and the 2dFGRS survey. Therefore, even though we 
use bias as a free fitting parameter, it is strongly constrained by the fact that 
the CMB and 2dFGRS data essentially cover much of the same range in /c-space. 

It should be noted that Elgar0y and Lahav find that bias does not play 
any role in determining the bound on '^m^. At first this seems to contradict 
the discussion here, and also what was found from a Fisher matrix analysis in 
Ref. 122 ■ The reason is that in Ref. [311, redshift distortions are included in 
the 2dFGRS data analysis. Given a constraint on the amplitude of fluctuations 
from CMB data, and a constraint on fimh^ , this effectively constrains the bias 
parameter. Therefore adding a further constraint on bias in their analysis does 
not change the results. 

Neutrinoless double beta decay 

Recently it was claimed that the Heidelberg-Moscow experiment yields pos- 
itive evidence for neutrinoless double beta decay. Such experiments probe the 
'effective electron neutrino mass rriee — \ ^ej'^i'j I- Given the uncertainties in 
the involved nuclear matrix elements the Heidelberg-Moscow result leads to a 
mass of rriee = 0.3 — 1.4 eV. If this is true then the mass eigenstates are necessar- 
ily degenerate, and X] "^j^ — Swee. Taking the WMAP result of ^ < 0.70 eV 
at face value seems to be inconsistent with the Heidelberg-Moscow result |47) . 
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However, already if hy-a forest data and a constraint on the bias parameter is 
not used in the analysis the upper bound of ^rrii, < 1.01 eV is still consistent. 
For this reason it is probably premature to rule out the claimed evidence for 
neutrinoless double beta decay. 

3.2 Neutrino relativistic energy density 

For the case of the effective number of neutrino species we have in all cases 
calculated constraints in the {flbh? , N^) plane, while marginalizing over all other 
parameters. The reason for this is that for Big Bang Nucleosynthesis purposes 
these are exactly the important parameters. Therefore, to combine CMB, LSS, 
and BBN constraints the marginalization over [2i,h^ should not be performed. 
Furthermore, when constraining N,^ we have always assumed that ^ ~ so 
that the neutrino mass has no influence on cosmology. 

We start out by investigating constraints on TV^ from CMB and LSS data 
alone. In Fig. 2 wc show Ax^ for a global fit to which marginalizes over all 
other cosmological parameters. The overall best fit for the WMAP T and TE 
data, combined with the Wang et al. compilation, the 2dFGRS data, the HST 
key project data on i/o, as well as the SNI-a data on i?™, has = 1467.6 for 
a total of 1395 degrees of freedom. This gives x^/d.o.f. = 1.05 which is entirely 
compatible with the best fit WMAP value for the standard ylCDM model of 
X^/d.o.f. = 1.066. We also show constraints for two other analyses. The first is 
for WMAP and 2dFGRS data alone and the second for WMAP data alone. The 
bounds for the three cases are 

= A.Qt2A for WMAP+Wang+2dFGRS+HST+SN-Ia 

= i.ltil for WMAP+2dFGRS 

= 2.ltil for WMAP only 

These bound are entirely compatible with those found by Grotty, Lesgourgues 
and Pastor jSOj, and much tighter than the pre- WMAP constraints. 

The constraints derived here are also compatible with what is found by Pier- 
paoli |3n|j where are assumption of spatial flatness was relaxed. 

In the lower panels of Fig. 2 we show the best fit values of i/o and i7,„ for a 
given value of N^. The main point to note is that the constraint on is strongly 
dependent on iJo- This was also found in Ref. With only GMB data and 
the weak top-hat prior on the bound on Ni, is very weak. Adding the HST 
Key Project prior on Hq cut away a significant amount of parameter space both 
at low and high N^. Adding the 2dFGRS and Wang et al. data mainly has the 
effect of shifting the best fit value to higher iVj^, but also cuts away the low N^, 
values, an effect also seen in Ref. |5U| . 

In Fig. 3 we show constraints on {fibh"^ ,N^) for the full data set described 
above. The best fit value for Qi,h^ is 0.0233, which is equivalent to the value 
found in the WMAP data analysis. In the 2-dimensional plots the 68% and 95% 
regions are formally defined by Ax^ = 2.30 and 6.17 respectively. Note that this 
means that the 68% and 95% contours are not necessarily equivalent to the same 
confidence level for single parameter estimates. 
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It should be noted here that in addition to an upper bound on Ni, there 
is also a S.Oct confidence detection of N^, > 0. This is in concordance with the 
pre-WMAP data from which a non-trivial lower bound on N^, could also be 
derived. 

Adding BBN information - In the case where all the relativistic energy 
density contained in Ni, is produced prior to BBN, a BBN constraint can be 
added to the CMB and LSS constraint without any problems. In practice we 
have used abundances of He-4 and D to make constraints in the (l?b/i^,7V^) 
plane. We use the following values for the primordial abundances |5SI54| D/H = 
2.78;g;3g X 10-5 and Yp = 0.238 ± 0.005 

This calculation is shown in Fig. 4. In terms of a single parameter constraint 
on Ni, it is iV^ = 2.6!to'3 (^5 % conf.). Compared to the recent calculation 
by Abazajian |S2 of a BBN-only constraint of 1.7 < N^, < 3.5 (95 % conf.) 
this is a significant improvement. Very interestingly the new limit suggests the 
possibility that N,y is actually less than 3. This is for instance possible in scenarios 
with extremely low reheating temperature |55I5()| . Note that this result is also 
consistent with the calculation by Barger et al. although t heir preferred 
region is slightly different because they do not include large scale structure data 
in their analysis. 

Of course this conclusion is mainly based on the fact that CMB and LSS data 
prefers a slightly higher value of Dbh^ than BBN. It should also be stressed that 
the estimates of the primordial abundances could be biased by systematic effects 
so that the quoted statistical error bar is not really meaningful. Therefore it is 
probably premature to talk of any inconsistency between the — 3 prediction 
of the standard model and observations. 

In fact the argument can also be reversed. If is fixed to the standard model 
value of 3 then then CMB and LSS constraint on flbh^ provides an accurate 
measurement of primordial He-4. Using the derived constraint on flbh? the 95% 
confidence range for Yp is 0.2458 <Yp< 0.2471. This could point to a serious 
underlying systematic effect in observational determinations of Yp, as discussed 
in Ref. |SH|- 

3.3 Combining ^ with AT^ - constraining LSND 

From the analyses in the above two sections it was found that: (a) An increasing 
^ nil, can be compensated by a decreasing and (h ) An increasing Ny can be 
compensated by an increasing Hq. One might therefore wonder whether a model 
with non-zero ^ , combined with Ni, > 3 can provide a good fit to the data. 

In order to test this we have performed a likelihood analysis for ^ rrii, for 
different values of N^. We show this in Fig. 5. This analysis was performed with 
all available data and priors. 

As can be seen from the figure, the best fit actually is actually shifted to 
higher ^ rrii, when N^, increases, and the conclusion is that a model with high 
neutrino mass and additional relativistic energy density can provide acceptable 
fits to the data. As a function of Ni, the upper bound on J2 is (at 95% 
confidence) 
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J2m^< 1-01 eV 
J2m^< 1.38 eV 
J2m^< 2.12 eV 



for N^^3 
for N^ = 4, 
for N^^5 



This has significant implications for attempts to constrain the LSND experi- 
ment using the present cosmofogical data. Pierce and Murayama conclude from 
the present MAP limit that the LSND result is excluded [H] (see also Ref. |1H]). 

However, for several reasons this conclusion does not follow trivially from 
the present data. In general the three mass differences implied by Solar, at- 
mospheric and the LSND neutrino measurements can be arranged into either 
2-1-2 or 3-1-1 schemes. Recent analyses HHj of experimental data have shown that 
the 2-1-2 models are ruled out. The 3-1-1 scheme with a single massive state, 
7714, which makes up the LSND mass gap, is still marginally allowed in a few 
small windows in the {Am? , sin^ 29) plane. These gaps are at {Am^, sin^ 29) ~ 



(0.8eV^2xl0-3),(l.8eV^8xl0-'*),(6eV^1.5xl0-3) and (10eV^ l.SxlQ-^). 



These four windows corresponds to masses of 0.9, 1.4, 2.5 and 3.2 eV respectively. 
From the Solar and atmospheric neutrino results the three light mass eigenstates 
contribute only about 0.1 eV of mass if they are hierarchical. This means that 
the sum of all mass eigenstate is close to 7714. 

The limit for N^, = 4 which corresponds roughly to the LSND scenario is 
^ TTii, < 1.4 eV, which still leaves the lowest of the remaining windows. The 
second window at to ~ 1.8 eV is disfavoured by the data, but not at very high 
significance. 

4 Discussion 

We have calculated improved constraints on neutrino masses and the cosmolog- 
ical relativistic energy density, using the new WMAP data together with data 
from the 2dFGRS galaxy survey. 

Using CMB and LSS data together with a prior from the HST key project on 
Hq yielded an upper bound of ^ rrii, < 1.01 eV (95% conf.). While this excludes 
most of the parameter range suggested by the claimed evidence for neutrinoless 
double beta decay in the Heidelberg-Moscow experiment, it seems premature to 
rule out this claim based on cosmological observations. 

Another issue where the cosmological upper bound on neutrino masses is very 
important is for the prospects of directly measuring neutrino masses in tritium 
endpoint measurements. The successor to the Mainz experiment, KATRIN, is 
designed to measure an electron neutrino mass of roughly 0.2 eV, or in terms of 
the sum of neutrino mass eigenstates, "^nii, < 0.75 eV. The WMAP result of 
^ < 0.7 eV (95% conf.) already seems to exclude a positive measurement of 
mass in KATRIN. However, this very tight limit depends on priors, as well as 
Ly-Qf forest data, and the more conservative present limit of '^rrii, < 1.01 eV 
(95% conf.) does not exclude that KATRIN will detect a neutrino mass. 

From the data we also inferred a limit on N^y of N^, = 4.012 ° (95% conf.) on 
the equivalent number of neutrino species. This is a marked improvement over 
the previous best limit of roughly N^, < 13 |51I60| . 
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When light element measurements of He-4 and D are included the bound 
is strengthened considerably to Njj — 2.6tQ 3 (^^ ^ conf.). Interestingly this 
suggests a possible value of which is less than 3. This could be the case 
for instance in scenarios with very low reheating temperature where neutrinos 
were never fully equilibrated [55*56'. However, it should be stressed that primor- 
dial abundances could be dominated by systematics. Therefore it is probably 
premature to talk of a new BBN "crisis" . 

Finally, we also found that the neutrino mass bound depends on the total 
number of light neutrino species. In scenarios with sterile neutrinos this is an 
important factor. For instance in 3+1 models the mass bound increases from 1.0 
eV to 1.4 eV, meaning that the LSND result is not ruled out by cosmological 
observations yet. 
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Fig. 1. The top panel shows as a function of ^ rriu for different choices of priors. The 
dotted line is for WMAP + 2dFGRS data alone, the dashed line is with the additional 
Wang et al. data. The full line is for additional HST and SNI-a priors as discussed 
in the text. The horizontal lines show Ax^ = 1 and 4 respectively. The middle panel 
shows the best fit values of Hq for a given ^m^. The horizontal lines show the HST 
key project la limit of Ho = 72±8 kmsMpc"^. Finally, the lower panel shows best fit 
values of fim- In this case the horizontal line corresponds to the SNI-a la upper limit 
of < 0.42. 




Fig. 2. as a function of ^ A^'^ for different clioices of priors. The dotted line is for 
WMAP data alone, the dashed line is with the additional Wang et al. and 2dFGRS 
data. The full line is for additional HST and SNI-a priors as discussed in the text. The 
horizontal lines show Ax^ = 1 and 4 respectively. The middle panel shows the best fit 
values of Hq for a given Ni,. The horizontal lines show the HST key project la limit of 
Ho = 72 ± 8 kmsMpc"^. Finally, the lower panel shows best fit values of Om- In this 
case the horizontal line corresponds to the SNI-a la upper limit of f2m < 0.42. 




Fig. 3. 68% and 95% confidence contours in tiie {^th^, N^) plane for tlie WMAP TT 
and TE data, combined witfi tfie 2dFGRS data, tlie HST data on Ho and tiie SNI-a 
data on Orn- 
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Fig. 4. 68% and 95% confidence contours in the {fifjh^ , N^) plane for the same data 
sets as in fig. 3, but with the addition ol BBN data. The lined contours are the 68% 
and 95% regions for BBN data alone. 
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Fig. 5. Ax^ as a function of ^ m„ for various different values of Ni,. The full line is 
for N„ = 3, the dotted for N„ = 4, and the dashed for A^i^ — 5. Ax^ is calculated 
relative to the best fit Nv = 3 model. 



